HANKEL DETERMINANT AND ORTHOGONAL POLYNOMIALS FOR 
THE GAUSSIAN WEIGHT WITH A JUMP 

A. ITS AND I. KRASOVSKY 

Dedicated to Percy Deift on the occasion of his 60th birthday. 

Abstract. We obtain asymptotics in n for the n-dimensional Hankel determinant whose 
symbol is the Gaussian multiplied by a step-like function. We use Riemann-Hilbert analysis 
of the related system of orthogonal polynomials to obtain our results. 



1. Introduction 
Consider the n-dimensional Hankel determinant 

(POO \ n — 1 POO POO n 

/ x j+k w(x)dx) =-/ •••/ Hixi-Xj) 2 Y[w(x k )dx k 
J—oc J j t k=0 n ' J-oc J -oo i< - k=1 

with discontinuous symbol 

(1.2) w(x) = e~* 2 S^e(-l/2,l/2). 

I e zpn , x > fi 

We are interested in the asymptotics of D n {j3) for large n. Note that if the jump in w(x) is 
absent, i.e. (3 = 0, D n (0) is a Selberg integral with the following explicit representation and 
the asymptotics (cf. [34]): 

n-l 

(1.3) D n (0) = (2it) n ' 2 2- n2 / 2 J] j\ = (27i) n (n/2) n2 / 2 n- 1 / 12 e- (3 ^ n2+ <'(- 1 \l + 0(l/n)), 

3=1 

where C(x) is the derivative of Riemann's zeta-function. 

The jump discontinuity similar to the one in (jl.2p and root-type singularities of the form 
\x — fi\ a are collectively known as Fisher-Hartwig singularities. 

Investigation of asymptotic behaviour of general Toeplitz determinants whose symbol has 
such singularities was initiated by Lenard [29] and Fisher and Hartwig |19j . who conjectured 
asymptotic formulas on the basis of explicitly known examples and Szego's theorem for 
smooth symbols. (Such determinants are related, for example, to the Ising model and to 
random walks on a lattice.) A proof of part of the Fisher-Hartwig conjectures was obtained 
by Widom in 1973 [35]. Since then many workers have contributed to proof /disproof of these 
conjectures (see [17] for a review) using mostly methods from operator theory. 

In contrast, the question of asymptotics for Hankel determinants, especially when the 
symbol has unbounded support, remains to a large extent open. Studies for in-a-sense 
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regular perturbations of Hankel symbols were carried out in [23| El EHl H] . The investigation 
of the "singular" Fisher-Hartwig case for symbols on R started with the particular situation 
of the Gaussian perturbed by root-type singularities as this example is important for random 
matrix theory; namely, for symbols 



1 2a, 



(1.4) e-* 2 H\x- H \ 

3=1 

It turns out that the proper singular case corresponds to fij being inside the support 
of the equilibrium measure, namely, [ij G [— \/2n, V2n). The asymptotics of the Hankel 
determinant with the symbol (11. 4p and fij = \j\/2n, Xj G (—1,1) were found in 0, ETJ [22] 
for integer <x, , and in for the general case 9ft<x, > —1/2. 

In the present paper we consider another (jump) type of Fisher-Hartwig singularity on the 
Gaussian background. We prove 



Theorem 1.1. Fix A G (-1, 1), 3?/3 G (-1/4, 1/4). Let fi = X V2n. Then, as n -> oo, 

= G(l + 0)G{1 Xlr^ /2 (8n)^ x 
(1.5) n[ J , , . , r 

'Inn 



1 + 



1-4|J?/3| 



n 



exp ^2in/3 ^arcsin Ao + Ao-y/l — Aq^ 

where G(z) is Barnes' G-function. 

With increasing effort, one can compute higher order asymptotics. The restriction $1(3 G 
(—1/4, 1/4) only makes it easier to handle some technicalities. One could use our methods 
to consider the general case G (—1/2, 1/2] as well. 

Note that 

(1.6) sp + W-«=(g)'«p(-^-f^), 

where T(z) is Euler's T-function. 

It is interesting to compare Theorem 11.11 with the result of [21 [8] on a Toeplitz determinant 
for a symbol with a jump. We have the same combination of G-functions in both cases. 

Of independent interest is the investigation of the system of polynomials orthonormal 
w.r.t. the weight ( 11. 21) . Consider the system Pk(x) = HkX k + ■ ■ • , ^ 0, k = 0, 1, ... , 
satisfying 

POO 

(1.7) / pj(x)p k (x)w(x)dx = 6jk, k = 0, 1, 



5 

— OO 



This orthogonality relation is equivalent to 

/°° 5 
Pk{x)x :j w{x)dx = — , j — 0, 1, . . . , k, k = 0, 1, . . . 
-oo H k 

We also consider monic polynomials pk(x) = pk(x)/xk = x k + ■ ■ • , which can be defined by 
the conditions 



1.9) / p k (x)x J w(x)dx = h k 6 jk , h k ^0, j = 0, 1, . . . , k, k = 0, 1, . . . 
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Obviously, x^T 2 = h^. 

The existence of such a system of orthogonal polynomials with nonzero normalization 
coefficients hk is well known for real weights. If 5ft/? 7^ 0, we will show existence of these 
polynomials for large enough degree. 

As any system of orthogonal polynomials, p n (x) satisfy a three-term recurrence relation: 



;i.io) 



Xp n (x) = p n +l{z) + Anp n (x) + B n p n _ x {x), 



n = 0, 1, . . . 



P-i 



0. 



Multiplying (11.101) by p n _i(x)w(x) and integrating gives 



'1.11^ 



B„ 



K-l / Pn(x)p n -l(x)xw(x)dx 



x; 



n-l 



K-l X 2 n 

Let us fix the notation for the two leading coefficients of the polynomials p n (x) as follows: 

Pn (x) =x n + p n x n ~ l + ln x n ~ 2 + ■■■ 



Comparing the coefficients at the powers x n and x n 1 in the recurrence relation (11.101) . we 
obtain the following identities we use later on: 



(1.12) 



.4, 



Pn Pn+1: 



X, 



n-l 



X, 



In - 7n+l - Pn + PnPn+1- 



Our methods allow us to compute various asymptotics for polynomials p n {x) as n — > 00. 
One such question is how the singularity of the weight affects the recurrence coefficients 
A n , B n for large n. This was discussed by Chen and Pruessner [TU] who found a behavior 
oscillating in n involving cos(cilnn + C2), but did not determine the constant phase Ci- We 
prove 



Theorem 1.2. Fix A G (-1,1), G (-1/2,1/2), p ^ 0. Let fi = X V2n. Then as 
n — > 00 ; the recurrence coefficients satisfy 



A, 



ip 2 sin Ti P 



;i.i3) 



((a + a- 1 )(8n) /3 (l - A 2 ) 3/3/2 r(-/3)e^+ (Ao) 



47rV2n 

+i(a - a- 1 )(8n)- /3 (l - \ 2 )- 3 ^ 2 T(p)e- nMXo) y 



O 



1 



n l-2|5R/3| 



;i.i4) 



where 



B„ = — 



ip\ c 



7rsinh0 + (Ac 



2 2 V /T^A 2 4(l-A 2 )sin7r/5 
x [(8n) 2/3 (l - A 2 ) 3 ' 3 r(/3)- 2 e (2n+1) * +(Ao)+iarcsinAo 

+ (8n)- 2/3 (l - A 2 )- 3 ' 3 r(-/3)- 2 e- (2n+1 ^+ (Ao) - iarcsinAo ] + O I — 



1 



n l-2|SR/3| 



a(z) 



1- Aq 

1 + Ao 



1/4 



.,i7r/4 



0+(Ao) 



2z / VT~— x 2 <ir. 
A 
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In particular, for (3 = ij, 7 6 R \ {0} ; 

A n = - =2 = sin 2 (n (arcsinAo + A \/l - A§ - ^ ] - 7ln(8n) 



;i.i5) 



;i.i6) 



V2n(l - A 2 ) 

+ arg r(z 7 ) - ^ ln(l - A 2 ) + arctan ^== | I 0( 1 / „ i | ^ 

2 1 + Vl-A^ 

B n = - + ^== + ^— ^ cos ^arcsinAo + Ao^l - A 2 

x cos ^(2n + 1) ^arcsin A + A ^/l — A 2 , — ^ — 27 ln(8n) 
+2argr(i 7 ) -37ln(l-Ag) - arcsinAo) +0{l/n). 



If furthermore Xq = 0, the above equations become 

;i.l7) A n = sin 2 (7ln(8n) - ar g r(i 7 ) + im/2)[l + 0(1/ n)}, 

V 2n 



(1.18) B n = - + 1 cos(2 7 ln(8n) - 2 ar g r(i 7 ) + vr(n + 1/2)) + 0(l/n). 

Except for the phase argr(ry), (11.171) and (jl.lBp were found in [TO]. Note also that similar 
formulas for the recurrence coefficients for the Jacobi weight perturbed by a jump were 
conjectured by Magnus [SO]. Our approach can be modified to verify this conjectureEl We 
will address this question in a future publication. 

Our proofs of Theorems 11.11 and 11.21 are based on the Riemann-Hilbert problem (RHP) 
approach to orthogonal polynomials [20] (see the next section), combined with the steepest 
descent method for RHP of Deift and Zhou [13] which proved to be very successful in the 
analysis of asymptotic behavior of orthogonal polynomials and related determinants (see, 
[TT] for an introduction and bibliography of earlier works, and [2S1 [2Z1 E31 ESI El E21 EH EH 
1251 d51 EE])- The steepest descent method for a RHP for the orthogonal polynomials given 
by expression (jl.7p allows us to prove Theorem ll.2[ The main technical difficulty here is the 
construction of an approximate solution (parametrix) of a related RHP in the neighborhood 
of An. We find a new parametrix given in terms of the confluent hypergeometric function. 

In order to prove Theorem 11.11 an additional analysis is needed. It is based on a classical 
formula connecting D n ((3) and the orthogonal polynomials Pk(z) = XkZ k + ■ • ■ , namely 



(1.19) D n (J3) = l[. 

j=o 

From this expression (cf. [HI O [2U [23 [16]) one can derive an identity for (d/d(3) hiD n (f3) 
in terms of only p n (Xo), Pn-i(Ao), and several leading coefficients of p n (x). Substituting the 
asymptotic expansions for the polynomials in this identity and integrating over (3 from (3 = 
to some (3, we obtain the asymptotics of D n (0)/D n (O). 



1 In fact, Magnus conjectured asymptotics of the recurrence coefficients for a perturbation which is a jump 
multiplied by a root-like singularity at the same point. Our approach would work in this general case as well. 
The case of the Jacobi weight perturbed by only root-like singularities was resolved by Vanlessen in [33j . 
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2. RlEMANN-HlLBERT PROBLEM FOR Pk(z) 

Consider the following Riemann-Hilbert problem for a 2 x 2 matrix valued function 
Y(z, k) = Y(z) and the weight w(x) given in (II. 2ft : 

(a) Y(z) is analytic for zGC\f. 

(b) Let iel\ fi . Y has L 2 (1R) boundary values Y + (x) as z approaches x from above, 
and Y-(x), from below, related by the jump condition 

(2.1) Y + (x) = Y_(x)(l "VP), xeR\» . 

(c) Y{z) has the following asymptotic behavior at infinity: 

(2.2) Y(z) = (l + o(-)) ( Z * z ° k ) , as 2 - oo. 



It is easy to verify that, if the polynomials satisfying (11.91) exist, this problem has a solution 
given by the function: 



(2.3) Y(z) 

-OO f- 

Note, in particular, that detY(,z:) = 1. Indeed, from the conditions on Y(z), detY(z) is 
analytic across the real axis, has all singularities removable, and tends to 1 as z — > oo. It 
is then identically 1 by Liouville theorem. Also, the solution is unique: if there is another 
solution Y(z), we easily obtain by Liouville theorem that Y(z)Y(z)^ 1 = 1. 

In the next section we derive an expression for (d/dj3) In D n (f3) in terms of the matrix 
elements of (12.31) . and in the section after that we compute the asymptotics of Y(z) using 
(a) -(c). 

The existence of the system of orthogonal polynomials Pk(z) = x k z k + ■ • ■ satisfying 
(ll.7p with nonzero leading coefficients x k for real w(x) is a classical fact. Moreover, the 
coefficients = D k /D k+ \ are regular functions of (3 (as follows from the determinantal 
representation for D k [32], E]). For all complex (3 in any fixed closed bounded set of the 
strip —1/2 < \Rf3 < 1/2, (denote this set by Q), we shall prove below the existence of a 
solution to the Riemann-Hilbert problem for all n larger then some n > 0. Its asymptotics 
will be explicitly constructed. We shall also see that the coefficients x k , k > n Q are nonzero 
and finite for all such (3. For k < no, the coefficients h\ are regular and nonzero (as follows 
from the determinantal representation) outside of a possible subset Q of fl. As a consequence 
of the determinantal representation, the system of the orthogonal polynomials exists, and 
the formula (11.191) holds for (3 G f2 \ Q. Throughout Section 3 (and, hence, in the differential 
identity obtained there) we assume that f3 Q. (There is no such condition on (3 in Section 
4.) This provides in Section 5 a proof of Theorem II .11 for (3 outside the set fl. However, as we 
shall see below, the error term in the asymptotics of D n is uniform for all (3 £ fl. Theorem 
11.11 will follow then by continuity. 

3. Differential identity 

We now derive a differential identity needed to prove Theorem 11.11 Throughout this 
section, we consider n a fixed positive integer and (3 G fl \ fl (see previous section). We start 
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with the following identity obtained in ([25], eq. (17)) for a weight w(x) depending on a 
parameter (3: 



J i = / V f n ^(x)i/ n _ lia! (x)w(x)dx, J 2 = p' ntX (x)p' n _ l0 (x)w(x)dx. 

J— oo J —oo 

Here the prime and the lower index f3 (respectively, x) stand for the derivative w.r.t. (3 
(respectively, x). 

This identity contains polynomials of order n and n — 1 only and therefore can be used 
to calculate asymptotics of Hankel determinants for rather general weights w(x). However, 
as in the situation for weights with root-type singularities considered in [25], the weights 
with jumps, as will be shown below, allow us to reduce the above identity to a "local" 
one involving only the values of the matrix elements of Y(z,n) at particular points. More 
precisely, for our weight (11.21) we will obtain an identity in terms of p n (^o), p n -i(p*o), their 
derivatives, and several polynomial coefficients. 

Note that p n (x) are analytic functions of (3, (3 Q, as follows, e.g., from their representa- 
tion as a determinant. 

For the reader's convenience, we recall first the derivation of (13. ip . Using (11.19P and the or- 
thogonality condition (11.81) . we have (note that Pj^ix) = ^j^x 3 ^polynomial of degree less than j) 



-lnD n (/3) = ^InJI^ 2 = " 2 E^ = " 2 E / Pi(x)^(x)w(x)dx 

(3.2) 



/]pj( x ) I w{x)dx. 



By the Christoffel-Darboux formula (e.g., 



n-l 

(3.3) y2p% x ) = — (Pn^^K-lW - Pn(x)Pn-l,xi X ))- 



3=0 



Substituting (13. 3p into (13. 2p . differentiating it w.r.t. /3, and using orthogonality, we obtain 
Let us now evaluate J\. Integrating by parts, 



Jl= PnA X )Pn-l,xi X ) e X e inf3 dx+ / PnA X )Pn-lA X ) e * e ^ 

J —OO J jjLQ 

_ 2 f°° 

(3.4) = 2sinh(z7r/3)e ^'^^(/ioK-^/io) - / p n -i(x)p' nAx (x)w(x)dx 

J -co 

POO 

+2 / p n ^i(x)p n j 3 (x)xw(x)dx. 



X 
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The first integral in the r.h.s. equals nx^ /3 /x n _i by orthogonality, and the second can be 
written as 



B 



p n _i(x)(x^ )/3 x rt+1 + (x n p n )'px n + (x n 7„)^r n 1 H )w(x)dx. 



As in (25], expanding x n+l and x n in terms of the polynomials Pj(x), we obtain, by orthog- 
onality and (11.121) . 



(3.5) 

Finally, 
(3.6) 



B 



x„ 



,y / \ 2 

K n,f3 ( Xn-1 



X„ 



X„ 



J\ = —n- 



x' 



n,/3 



+ 2sinh(i7r/3)e ^^(/ioK-iM + 25. 



A similar calculation for J 2 yields 

x 7 

j 2 = 2 ^M + 2sinh(m/5)e-^ n (/i K_ li/3 (/i ). 



(3.7) 



x, 



(3.8) 



Substituting the above expressions into (13.11) and using (12.31) . we obtain 

Proposition 3.1. Let (5, —1/2 < 3?/5 < 1/2, be such that the system of orthogonal polyno- 
mials Pk{z) with finite non-vanishing leading coefficients satisfying ( [1. 7\ ) exists. Fix n > 1. 
Let p n = x n (x n + (3 n x n ~ 1 + ■y n x n ~ 2 + •••), and the matrix Y(z) be given by Ii2.3\) . Then 

d ( x 2 \ ' 

— In D n ({3) = -n(ln x n x„_i)^ - + 2 (7^ - f3 n {3' nt(3 ) + 

1 2 

-e~"° sin (7T/3) [Y n (n )Y 2 i(no)p - K u (// ) ^21(^0) - (In ^Xn-i)J*ii(//o)*2i (/*<))] • 

4. Asymptotic analysis of the Riemann-Hilbert problem 

The standard steps of the steepest descent analysis applied to the Riemann-Hilbert prob- 
lem (a) - (c) of Section 2 involve rescaling by \/2n so that the support of the equilibrium 
measure becomes [—1,1], deformation of the contour such that the jump matrix becomes close 
to the identity, construction of local parametrices around the points 1,-1, and Ao = /io/ V2n, 
and matching the parametrices with the solution in the region outside those points. All these 
steps are standard except for the construction of the parametrix around Ao- In this section 
we carry out the analysis and obtain Theorem 11.21 as well as the asymptotics necessary for 
the proof of Theorem 11.11 in the last section. Henceforth we assume that Ao G (—1, 1). 

4.1. Three transformations of the Riemann-Hilbert problem. We now perform stan- 
dard transformations of the RHP. The first one Y — > U is a scaling: 



(4.1) 



Y{zV2^) = {2n) na3/2 U{z), 



0-3 



1 
-1 



The second one U — > T is given by the formula 



(4.2) 



U{ Z ) = e ^/2 T ^ e n(g(z)-l/2)a^ 



cS 
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where 

I = -1 - 2 In 2, 



(4.3) 



g(z) — / hi(z — s)ip(s)ds, z & C \ (— oo, 1], ?/>(z) = —VI — z 2 . 
J-i * 

Below we always take the principal branch of the logarithm and roots. The function ip(z) is 
the scaled asymptotic density of zeros of p n (z). (Note that J^ 1 ^jj(x)dx = 1.) The function 
g(z) has the following useful properties: 

g + (x) +g-(x) - 2x 2 - I = 0, for x e (—1,1) 
g+(x) +g-(x) -2x 2 -I < 0, for igM\[-1,1] 
(4-4) ( 2ni, for x < -1 

9+(x) ~ 9-(x) = I 2vri £ ijj(y)dy, for x G [-1, 1] 
I 0, for x > 1 

From the RHP for Y(z), we obtain the following problem for T(z): 

(a) T(z) is analytic for z G C \ R. 

(b) The boundary values of T(z) are related by the jump condition 

rp f ■ s _r r f s / e -"to+(*)-ff-(*)) w (z) \ , 1 n 

J + l x J — 1 -\ x ) I q e n( 9+ (a:)-g_(a;)) J > x fc l — i J- 

( L )) t f \ r ( s fl u(x)e n ^)+9-i-)-^-i)\ 

T +(x) =T-(x) I w 1 I, a; G R\ [-1,1]. 

(c) T(z) = J + 0(l/z) as z -> oo, 

where are the values on the real axis of the following function defined in the whole 
complex plane with a slit along the line = A : 



(4.6) u(z) 



e lf3n , $iz < A 
e~ il3n , Kz > A ' 



Note that this problem is normalized to 1 at infinity, and the jump matrix on (— oo, — 1) U 
(1, oo) is exponentially close to the identity (see f)4.4p ). In order to have uniform bounds, we 
have to exclude small neighborhoods of the points —1 and 1, where g+(x) + g~(x) — 2x 2 — I 
is close to zero. 

Let h(z) be the analytic continuation of 

(4.7) h(x) = g+(x) - g_(x) = 2vre 3 ^ 2 f ^(y)dy 



to C \ ((— oo, — 1] U [l,oo)). A simple analysis shows that 9?/i(z) > for Sz > 0, and 
3?/t(z) < for 3z < in some neighborhood of (—1, 1). We again exclude neighborhoods of 
the points —1 and 1 for a uniform estimate. 

By the steepest descent method of Deift and Zhou, we now split the contour as shown in 
Figure 1. Define a new transformation of our matrix-valued function as follows: 
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Figure 1. Contour for the S'-Riemann-Hilbert problem. 



(4.8) S(z) = { 



(T(z), 
T(z) 



T(z) 



1 

-uj(z)~ 1 e~ nh(z ^ 1 

1 o N 



for z outside the lenses, 
, for z in the upper part of the lenses, 

for z in the lower part of the lenses. 



Then the Riemann-Hilbert problem for S(z) is the following: 

(a) S(z) is analytic for z E C \ E, where S = 1U LlJ =1 Iy. 

(b) The boundary values of S(z) are related by the jump condition 



S+(x) = S-(x) 







x e u. 4 =1 E„ 



(4.9) 



where the plus sign in the exponent is on E 2 , S 4 , and minus, on E 1; E 3 , 

u{xf 



S + (X) = S-(X) y 







xe (-1,1). 



S+(x) = S-(x) 



1 c<j(a;)e n ^ + ^^ )+!? "( :E )~ 2:E2 ~^ 
1 



x e M\ [-1,1]. 



(c) S{z) =1 + 0(1/ z) as z -> oo. 

Recalling the remarks above, we see that, outside the neighborhoods £/a > f±i, the jump 
matrix on Ej, j = 1, . . . , 4 is uniformly exponentially close to the identity. (Note that we 
had to contract the lenses at A , because of the cut of u>(z) across the line $tz = A .) We 
shall now construct the parametrices in U\ , U±i, and = C \ (U\ U XJ\ U C/_i). After 
that we shall match them on the boundaries dU\ and dU±i, which will yield the desired 
asymptotics. 

4.2. Parametrix in U^. We expect the following problem for the parametrix P^ in U^: 

(a) Poo (z) is analytic for z e C \ [— 1, 1], 

(b) with the jump condition on (—1, 1) 



(4.10) 



P (r)( ° 



xG (-1,1), 



(c) and the following behavior at infinity 



(4.11) 



Poo(z) — I + O ( - ) , as z ^ oo. 
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A solution Poo(z) can be found in the same way as, e.g., in [28] : 



(4.12) P M = i(D.r il^l-r 



-1 



-i[a — a 



-v 



V{z) 



-cr 3 



a[z\ 



z - 1 



a + a ± ) ~ v ~' \z + l 

where the cut of the root is the interval (—1,1). Note that det Poo{z) = 1. Here 



1/4 



(4.13) 



V(z) = exp 



^fz^l f 1 lno;(0 d£ 



2tt 



-i vi^e z -z 



Poo = lim £>(z) 



The Szego function T>(z) is analytic outside the interval [—1, 1] with boundary values 
satisfying V + (x)V_(x) = u(x), x G (—1, 1). Calculation of the elementary integral in (14.131) 
gives 



V(z) 



z\ -l-iy/(z2-l)(l-\l) 







(4.14) 



z - A 
z - A 



3 i7T/3/2 



^0-1 + V^ 2 - 1 )^-^). 
From here it is easy to obtain as the main term of the expansion in 1/z at infinity 



(4.15) 



&An + \/l — An ^ -^™A 



Below we will need 2 more terms in the expansion: 



(4.16) m 



(V(l - A|) + W„y/ 1 - Ajj) + 0(.-- 3 ). 



OO. 



At the point A , 

(4.17) ^)=c(/5)(z-A )' 3 (l + o(l)), c(/3) = 2^(1 - A^e"^ 2 , A , 3* > 0. 

Note that in order to obtain the correct argument of the complex number c(/3), we returned 
to the analysis of the original integral in (14.131) . 

Moreover, we will need the following expansions at ±1 which are easy to obtain from 
( I4.14p and the definition of a(z)- 



V\z) to(z 



+ 



uo(z) V 2 (z) 
V 2 (z) u(z) 



(4.18) u{z) V 2 (z) 



= 2(l -A^^u + 0(u 3 ^ 



a 2 + a 



a »- a -' = -Jl + y» + o<«n 



u = z — 1, 



z -> 1. 
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and 



V 2 (z) u(z) f „ 9 1-A nl o /9 A 
—FT- + = 21 + 4/3 2 - At + 0(u 3/2 ) 



^-#v = - 25/ Wl+r^ +0 <" 3/2 >- 

(4.19) u\ z ) v ( z ) Vl + Ao 

° 2 + a " 2 - > ~ Mf) + °(" s/2 »- ° 2 - «- 2 = 4 ()/! + 5 VI) + °(" 3/2 )- 

it = z + 1, Z — > — 1. 

4.3. Parametrix at the jump point. Let us now construct the parametrix P\ {z) in Z7> . 
We look for an analytic matrix- valued function in a neighborhood of U\ which satisfies the 
same jump conditions as S(z) on E n U\ , has the same behavior as z — > Ao, and satisfies 
the matching condition 

(4.20) P Xo {z)P-\z) = / + (1) 

uniformly on the boundary dU\ as n — > oo. 

Using the analytic continuation of V'(y) ( see (14. 3p ). define: 



(4.21) <f,(z) ~- 

Clearly, e^ z > is analytic outside [—1,1]. Set 



h(z)/2 = e 3in / 2 -Kf*ip(y)dy, for %z > 0, 
e-™h(z)/2 = e m/2 n j*i)(y)dy, for < ' 



(4.22) /(*) = 2vre^ 2 / ^{y)dy. 

Let us now choose the exact form of the cuts S in U\ j so that their images under the mapping 

C = n f( z ) are straight lines (Figure 2). Note that ((z) = nf(z) is analytic and one-to-one 
in the neighborhood of U\ , and it takes the real axis to the imaginary axis. We have 

(4.23) C = nf(z) = 4ne i7r / 2 ^l-X 2 (z - A )(l + 0(z - A )), z - A . 
We look for P\ (z) in the form 

(4.24) P Xo (z) = E{z)P {1 \z)e- n ^ z)a \ 

where E(z) is analytic and invertible in the neighborhood of U\ , and therefore does not 
affect the jump and analyticity conditions. It is chosen so that the matching condition 
is satisfied. It is easy to verify that P^\z) satisfies jump conditions with constant jump 
matrices. Set 

(4.25) P«(z) = * /J (C) = 

where ^ p{Q satisfies a RHP along cuts given in Figure 2: 
(a) m p is analytic for ( G C \ Uj^iy 
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III + 

r 4 

Figure 2. The auxiliary contour for the parametrix at A . 
(b) \Ijg satisfies the following jump conditions: 



(4.26) 




-(C) 


= *A-(C) 


(4.27) 




-(C) 


= *A-(C) 


(4.28) 




-(C) 


= *A-(C) 


(4.29) 


*/3,H 


-(C) 


= *A-(C) 







-17T/3 







for c e i\, 

_ e -i7r/3 o y ' for C e r 4 , 
for (er 2 u r 6 , 



l o 
l o 

g-iTT/3 X 



for c e r 3 u r 5 . 



We will solve this problem explicitely in terms of the confluent hypergeometric function, 
ip(a, c; z) for the parameters a = f3 and c = 1. For the reader's convenience we briefly sketch 
the standard theory of the function if) (a, c; z) in the appendix. 

Introducing the notation 

(4.30) ip(a,z) =ip(a,l;z), 

we define the following piecewise analytic matrix-valued function on the punctured complex 
plane C G C \ {0}: 

i/j((3, Qe 2mf5 e-V 2 -^(1 - 0, e-^Oe^e^ 2 ^ 

>(1 + (3, Oe^e-^M e"* C)e c/2 

< argC < 7r, 



(4.31) 



*(C) = 



for > 0, and 

*(C) 



W,C)e 



27ri/3 e -C/2 



(4.32) 



x 



o • • a i i 7r < arg C < 27r, 
^2? sin n p 1J & ^ ' 



^(l-/3,e— C)e™V/ 2 q^f 
V>(-/3, e~ i7r C)e c/2 
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for S£ < 0. @ Denote by Roman numerals the sectors among the cuts in Figure 2. We 
have the following 

Proposition 4.1. A solution to the above RHP (a), (b) for ^ is given by the following 
formulas: 

forCelU III, 
1 N 

-e inf3 " 

\ _ e -i?r/3£ N 
e m/3e q 

giTT/3 Q 



(4.33) 


*/j(C) 


= *(C) 


(4.34) 


MO 


= *(C) 


(4.35) 


*/*(C) 


= *(C) 


(4.36) 


*/j(C) 


= *(C) 


(4.37) 


*/j(C) 


= *(C) 



for ( e n, 

for C e IV, 
for CeV, 
for C e VI, 



where 



1, > o 
-l, 3C<o 



Proof. Equations (14.331) - (14.371) yield the correct jumps across all the contours Tj for any 
function \I/(C) which is continuous in the upper and lower ^-planes. Let us show that the 
choice (I4.3ip - (14.321) implies the continuity of ^(C) across the real axis. In the case of the 
negative semi-axes, this is trivial. Indeed, for ( < 0, in both (I4.3ip and (14.321) we have to 
take the same value of arg£, namely argz = ir, and hence we conclude at once that 

(4.38) 'MO = *-K>(-2isW ?)• for < <0 

(the real axis is oriented from the left to the right). This relation, in view of (I4.34p . yields 
the continuity of ^(C) across the negative semi-axes. 

It is in order to establish the continuity of ^(C) f° r C > 0, that we need the definition 
(I4.30p of ip((3, C) as a confluent hypergeometric function. The following two fundamental 
properties of the function ip{f3, C) will play the central role: 

27T? 

(4.39) m e 2m C) = e- 2CT V(/9, C) + e^™^(l ~ P> ^0^, 
and, by replacing (3 with 1 — (3 and e i7r £ with (, 

(4.40) ^(1 - 13, e™() = e 2i ^(l - /?, e^Q - e^— -^(/3, CK C . 

r 2 (i -p) 

These equations hold on the universal covering of the punctured £-plane, i.e. for all values of 
arg£, and they follow from the general relation (I7.30P proven, for the reader's convenience, 
in the appendix. 



2 We emphasize that, while the confluent hypergeometric function ip(f3, £) is defined on the universal 
covering of the punctured ( - plane (see the appendix), the function is defined on the punctured £ - 
plane itself. Thus the indication of the ranges of the argument of £ in the equations (|4.31|) and (|4.32p . is 
only relevant for the right-hand sides of the equations. 



14 A. ITS AND I. KRASOVSKY 

In the case ( > 0, the equation \E f ( g i+ (^) = is equivalent to the equation 

(4.4D * + (CH*-(o( c2 r - 2 :s^ 

(the real axis is oriented from the left to the right). For the entry (11) we have, 

^n_(C)e 2 ^ = W,C)e 4 ^e- f/2 

(4 42) TYl — 6) 

-2i sin(7r/3)V>(l - (3, e~ i7r C)e 37rt/3 e c/2 K H) , arg C = 2tt, 

J- IP J 

or 

^ n _(C)e 2i7r/3 = e 2 ^C)e 4m/3 e- f/2 
(4 ' 43) -2i sin(7r/3)V(l - /?, e l7r C)e 3m/3 e^ 2 ^^, arg C = 0. 

Eliminating, with the help of (14.391) . the function ip(/3, e 2m Q from the last equation, and 
taking into account the classical relation sin7rs = 7r/(r(s)r(l — s)), we obtain that 

^n-(C)e 2m/3 = MP, Oe 2 ^e- c/2 + e 3in(3 -^-^(1 - [3, e^C)e C/2 

(4 - 44) rri - T 

which, in view of the definition (14.31 j) . coincides with \I/ii + (£) for ( > 0. 
The entry (12) in the r.h.s. of (14.411) reads, 

-2zsin(7r/3)*n_(C) + ^-(CK 2 ^ 

(4.45) = -2i sin(7r/3)^(/3, C)e" C/2 - ^(1 - p\ e^K"^ 2 ^" ^ , arg C = 0. 

where we used expression (14.441) for ^ r n_(C), and the entry ^ 12 _(C) was taken from (14.321) 
with the substitution ( e 2m ( (as in the transition from (I4.42p to (14.431) ). 
Substituting (O01) for - /3, e in (), we obtain that (g3Sj> equals 

(4.46) - r( ^ ~ /?) e fa V/V(l - ft C), argC = 0, 

which is *&i2+(z) for ( > 0. 

Consider now the entry (21) of equation (14.411) . As for the entry (11), we first write the 
r.h.s. of the entry (21) as 

* 21 _( C )e 2 ^ = -^(1 + /?, e 2 -C)e 3 ^e~^ 2 ^^ 

+2i sin(7r/3)^(-/9, e^Qe^e^ 2 , arg C = 0. 
Now use (14.391) . with (3 replaced by 1 + /3 to eliminate ^(1 + /3, e 2ni ()'- 

* 2 i-(C)e 2 ^ = -e 3 ^e-^ 2 ^^ (V(l + ft C) 
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(4.48) 



(4.49) 



— e 



-17T/3 



2ni 



P»(l + /3) 



ip{-(3, e^C)e C ) + 2z sm(vr/?)^(-/5, e i7r C)e 27ri/3 e c/2 



. e 3Ti/3 e -C/2 



r(i + fl 
r(-/3) 



^(l + /3,C)e 



-2irift —iwf3 



2m 



2ni 



r 2 (i + /5) 

V>(-/3, e i7r C)e 2m/ V /2 



(4.50) 



, e ^ e -c/ 2 r(i +f) ^ (1 + AC); 



argC = 0, 



which, in view of the definition (14.311) . coincides with ^2i+{z) for £ > 0. 
Finally, for the entry (22) of f)4.4ip we have 



-2zsin(7r/3)* 21 _(C) + *22-(C)e 



(4.51) = 2i sin( 7 r/?)e-^e-^ r ] l|^V (l + ft C) + ^(-ft e^Oe" 2 ^ 2 , 



argC = 0, 



where we used expression ( 14.501) for ^21- (C)? an d the entry x E f 22— (C) was taken from ( 14.321) 
with the substitution ( -> e 2i7T C (cf- P5| ). Applying (r4~4"0"l) with /3 replaced by -(3 to 
P, e t7T (), we obtain that (14 . 5 1|) equals 



(4.52) e^(-/?,e-"C), argC = 0, 

which is \&22+ (z) for £ > 0. This completes the proof of Proposition 14.11 



□ 



We will now show that this solution can be matched with P^z) on the boundary dU\ 
for large n. We need first to compute the asymptotic expansion of ^(C) for C ~~ > 00 • F° r 
that we use the classical result [3] (see also the appendix) for the confluent hypergeometric 
function: 



(4.53) x) = x~ a - a l x~ a ~ Y + 0(x^ 2 ) 



\x\ 



00, — 37r/2 < argx < 37r/2. 



We can apply this result to (I4.31j) to obtain the asymptotics of the solution in sector /. 
The "correct" triangular structure of the right matrix factor in (14.341) implies that these 
asymptotics remain the same in the whole second quadrant, namely: 



(4.54) 



MC) = 

where we denoted 
(4.55) 



,27ri/3 







7T 



- < arg C < vr, 



/(*) 



r(i - o) 

T(a) 
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Similar considerations in sector VI and in the upper half of sector V, where we still can 
obtain the asymptotics by direct substitution of (I4.53P into (I4.3ip . yield 



(4.56) 



-/3 2 -f((3)e 



17T/3 



+ 0(C 



-2\ 



7T 



C -> oo, 0<argC<-. 



Before discussing the asymptotics of ^(C) i n the sectors in the lower half-plane, we need 
to perform the deck transformation, ( — > e _27 ™C, with the first column of the first matrix 
factor in the r.h.s. of (Q2)l . Using 



^(/3,C)e 2 " i/3 e- c/2 = ^(/3,e 27ri (e- 2wi C))e 2m/3 e 



-2l7T/3 



-17T/3 



2?ri 
f%3) 



^(l-/?,e- m C)eMe 27ri/3 e" c/2 



^(/?,e- 2 ^C)e" C/2 + e i7r/3 



2?ri 



^(l-/3,e- m C)e 



(4.57) 
and 



-^(1 + 



^ e- 2m C)e- ?/2 + 2* sin(vr/3)^(l - 0, e -- C )e^e^ 2 ^^ 
/3,C)e" /3 e-^ r ^ + f ) = -^(l + Ae 2 "(e- 2 "())e"V^ r(1 + /j) 



r(-/3) 

V(l + /?,e"^C)-e- 



r(-/3) 



2vn 



r 2 (i + /3) 



(4.58) 



-^(l + /5,e- 2m C)e~ m/3 e 



r(-/3) 

2rr '< w - ;T ' ; - ~ k/2 r ^ + ^ - 2isin(7r/?)^(-/5,e-^C)e C/2 . 

r( p) 



Equations (I4.57P and (14.581) allow us to rewrite (I4.32p in the form 

iP((3, e- 2 "C)e" c/2 - ft e~H)e^/ 2 ^=^ 



(4.59) *(C) 



-V(l + A e- 2 ^C)e-^e-^ 2 ^M 
7r < arg C < 27T. 



^(-/?,e" i7r C)e f/2 



Now the asymptotic result (14.531) can be applied directly for the whole range 7r < arg£ < 2tt, 
and as for (14.541) and (I4.56p . we arrive at the following asymptotics for the function ^(C) 
in the lower half-plane: 



(4.60) 



MO = 



I+ 1 - -^) e H+ofr 2 ) 



e- m P ) 



C -> oo, vr<argC<y, 
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and 



(4.61) 



*/j(C) = 

xC"^ 3 



-P 2 



-f(P)e 



iwf3 





-27TI/3 







^ 2 

C -> oo 



+ o(r 2 ) 

3tt 



< argC < 27T. 



Relations (I4.54p . (14.561) . (I4.60p . and (14.611) give a complete description of the asymptotic 
behavior of the function ^(C) m the neighborhood of ( = oo. In particular, for the left and 
the right half-planes we have that 
(4.62) 



*/»(C) 



/+-( 

71 37T 

- < argC < 



-f(P)e 
P 2 



in/3 



+ o(c 2 ) 



-/3<x 3 









—iirft 



-W2 



for §l( < half-plane (upper half-plane in z variable), and 



(4.63) 



<M0 = 



71 71 

~2 " argC - I" 



-/(/3)e 
/3 2 



17T/3 



D iri/3 



for -Ji^ > half-plane (lower half-plane in z variable). Substituting these asymptotics into 
the matching condition ( 14.201) : 

(4.64) Px (z)P-\z) = E(zMC)e- n ^p-\z) = I + o(l), 

we obtain 

' e -27Ti/3 Q 



(4.65) E{z) = P 00 (z)C Pa3 e nMXo)a3 

for Qz > 0, and 







^in(3 j i 



(4.66) 



E(z) = P 00 (z)C^ <T3 e 



/3o"3„-n0 + (A o )o- 3 



1 

_ e ~iKd Q 



7T 37T 

2 - argC -T' 



7T 7T 

"2 " argC - 2' 



for 3^ < 0. A direct check using the jump condition for Poo{z) shows that E(z) has no jump 
across the real axis. Moreover, E(z) has no singularity at z = Ao as easily follows from the 
expansion 

(4.67) ^ = (8n)^e-^(l - Aq)~ 3/3//2 (1 + o(l)), * - A , 3z > 0, 

which is derived from (14.171) . Therefore, is an analytic function in U\ . This completes 
construction of the parametrix at Ao: it is given by the formulas (14. 24)4. 25114. 65114. 66p and the 
formulas (14 . 33 p (14 . 3 T[) of Proposition 14.11 In particular, in the pre-image z(I) of the sector 
/ of the C-plane, we have: 

P Xn (z)=P 00 (z)(^x 



(4.68) 



^(/3, C) V(l - P, e - t7T C)f(p)e- inl3 e 2n ^ x ^ 



z e z(I), 
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where 

Y(l — a) f l , f z , 

(4.69) f(a) = \, . > , 0+(A o ) = 2i / y/l - x 2 dx, ( = Ain / Vl - x 2 dx. 

r 0) ' Jxo Jxo 

Note that the parametrix P\ (z) correctly reproduces the type of singularity of the solution 
of the RHP of Section 2 at Ao- Indeed, on the one hand, it is clear from (12. 3p that the solution 
has a logarithmic singularity at z = Xq (or ( = 0). On the other hand, the expansion of the 
function ip(a,l, () at zero is known to be [3] 

(4.70) if,{a, 1, C) = -pjU ( ln C + - 2^r) + O(Ch0, C - 0, 

where Cr = 0.5772. . . is Euler's constant. In fact, P\ {z) reproduces the singularity of the 
solution at A exactly, in the sense that S(z)P\ ( z ) 1 is analytic in a neighborhood of A . 
The latter fact is a consequence of the absence of jumps and the fact that the singularity 
can be at most logarithmic. 

We can extend (I4.20p into a full asymptotic series in n. For our calculations we need to 
know the first correction term: 

^Ao (z)P-\z) = I + AxCr) + 0(l/n 2 ~WI), 

(4.71) Ai(z) = ^Poo(z) L_2^ e -2n^+(Af)+27ri/9j^_^ 

z G dz(I), 

where dz(I) is the part of dU\ whose C-fmage is in /. As the calculation for the other sectors 
shows, this expression for Ai(^) extends by analytic continuation to the whole boundary dU\ 
(cf. [27] [33] [25]). Moreover, it gives rise to a meromorphic function in a neighborhood of 
U\ with a simple pole at z = A . The singularity given by actually cancels with that of 
T>(z) from P^ (see ( 14.671) ). The error term 0(l/n 2 ~ 2 ^^) in (14.711) is uniform in z on dU\ 
and in (3 in bounded sets of the strip — 1/2 < 3?/3 < 1/2. 



.^ e 2n0 + (A 3 -)-27ri/3 
P 2 



p: 



4.4. Parametrices at z = ±1. We now construct parametrices in the remaining regions 
JJ\ and U-i. These are obtained by a slight generalization of the results of [14J, which can 
be viewed as the case uj(z) = 1. The construction is identical with the corresponding one in 
[25]. 

We are looking for an analytic matrix-valued function in JJ\ which has the same jump 
relation as S(z) there and satisfies the matching condition on the boundary: 

(4.72) P 1 (z)p-\z) = I + 0(l/n). 

The solution is: 



P 1 = E(z)Q(0e- nv(z)as uj(z) 



(4 ' 73) E{z) = P 0O (^)cu(^) CT3/2 e^ 3/4 v / ^ (] A ^/4 e -Wi2 
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and Q(£) is given by the expression (7.9) of [H] in terms of Airy functions (in the notation 
of HU, Q(0 = ^ a (b)- In these formulas 

3 X2/3 



(4.74) £(z) = l-ne-^^z] 
is an analytic function in a neighborhood of z = 1 with a cut in the interval (—1,1), and 

(4.75) = 2n 2 /\z - 1) (l + 1(* - 1) + 0((z - l) 2 )) . 

The argument of the Airy function on dU\ is uniformly large, so we can expand it into an 
asymptotic series and proceed the same way as for dU\ . As a result we have the matching 
condition ( 14.721) extended to a full asymptotic expansion in the inverse powers of n. We shall 
need only the first 2 terms: 

P 1 (z)p- 1 (z)=I + A 1 (z)+0(l/n 2 ), 

(4.76) A 1 (z) = P^M^/V- / 4 1 _ 1 1 /6 ) e-^l^{z)-^p-\z)^ 
z e OUl 

The function Ai(z) is meromorphic in a neighborhood of Z7i with a second order pole at 
z = l. 

The argument for the parametrix in U-i is similar. We have: 

P_i = J E(z) ( T 3 Q(e^0^3e~ n ^ )<T3 w(2)- a3/2 , 

(4.77) . i*r„~iA I — / 1 1 



•3/2 



£7(2) = P 00 (^)cu(z) ff3/2 e i7r,J3/4 v / ^ f _^ J J (e -^ )CT3 /4 e — /i2 



Here 



(4.78) 



For z 



£(z) = e— ( -n0(z 



3 ~ . 



"/i(z)/2 — vk — e 3 ^/ 2 vr /* x ^(j/)dy, for 3z > 0, 
>™h(z)/2 + iix = e 5 ^/ 2 7r f* x ^(y)dy, for 3z < 



(4.79) = 2n 2 / 3 (l + (l - 1(1 + z) + 0((1 + z) 2 ) 

The first 2 terms in the matching condition: 

P- 1 (z)P oo 1 (z) = I + A 1 (z)+0(l/n 2 ), 

(4.80) = ^oo(^)^(^) CT3/2 e— / 4 1 ( X f _~J 6 ) e---/4 a;( ^ ) -.3/2p-i ( ^ ) 

As in (I4.7ip . the error terms in (14.761) (resp., (14.801) ) are uniform for all z G dU\ (resp. 
z G dU-x) and (3 in a bounded set. 
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U 




Figure 3. Contour for the i?-Riemann-Hirbert problem. 
4.5. Solving the RHP. Let 



(4.81) R(z) 



(S(z)P- 1 (z), 2£[/ 0O \S, 

S(z)P- o 1 (z), zeU Xo \Z, 

S(z)Pr\z), zeU!\E, 

[s(z)p:1(z), ze[f_i\s. 



It is easy to verify that this function has jumps only on dU±i, dU\ , and parts of 
M. \ [—1,1] lying outside the neighborhoods U±i, U\ (we denote these parts without the 
end-points E out ). The contour is shown in Figure 3. Outside of it, as a standard argument 
shows, R(z) is analytic. Note that R(z) = I + 0(1/ z) as z — > oo. 
The jumps are as follows: 

R+(x) = R.(x)Poo(x) ( L Q u[x)e 1 )P nci (x) 



(4.82) 



x eR\ [-1-5,1 + 5], 

R+(x) = R-(x)P 00 (x) ^^j-igTiftW ij P°°i x ) \ 



xGSr, fc = l,...,4, 

where the plus sign in the exponent is on E^f, and minus, on E^i-i, 3 = 1> 2> 
R + (x) = R-(x)P\ (x)P oo (x)~ , £ £ t^A \ { intersection points}, 
R+(x) = R-(x)P±i(x)P OQ (x)~ , x G <9£/±i \ { intersection points}. 

Here 5 is the radius of U\ and E7_i. The jump matrix on S out can be estimated uniformly 
in /3 as / + 0(exp(—en\x\)), where e is a positive constant. The jump matrices on dU\ 0i ±i 
admit a uniform expansion in the inverse powers of n multiplied by n 2 ' ^ ' (the first 2 terms 
of which are given by (H~7I1) . (pTTS]) . and (Mill ): 



(4.83) / + Ai(z) + A 2 (z) + • • • + A fc (z) + 0(n 



Every Aj is of order 7i~^+ 2 W. (The above expressions give us the explicit form of Ai(z) in 
each of the neighborhoods.) 

The solution to the above RHP is given by a standard analysis (see, e.g., [T4" t [15l [T6| 125]). 
We obtain that 

k 

(4.84) R{z) =I + J2 R ^ z ) + 0(n- fc - 1+2 l^l), i2,-(z) = 0(n^'+ 2|Ji/31 ), j = 1, 2 . . . 

3=1 
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uniformly for all z and for j3 in bounded sets of the strip — 1/2 < < 1/2. The expressions 
for Rj are computed recursively. We shall need only the first one: 

(4.85) RJ Z ) = -L [ Al(x)cb , dU = dU 1 UdU^ 1 UdU Xo . 

2m J du x - z 

The contours are traversed in the negative direction. 

4.6. Asymptotics for the polynomials and recurrence coefficients. The asymptotic 
solution given above allows us to obtain various explicit formulas for asymptotics of the 
polynomials orthogonal w.r.t. the weight ( 11. 2p . First of all let us note that such polynomials 
exist for large enough n, i.e., the normalizing coefficients hf. are nonzero. Indeed, tracing the 
transformations Y U ^ T — > S of the RHP, we obtain from (12.31) 

K\ = hm = hm r al ( g )e^W-0 

(4.86) w ; 

iV2n 

27i(2n) 



lim —^—zSv(z)e- nl . 



z— >oo 



Note that S(z) = (J + 0(l/n l 2 ^^))P (X (z); and a simple calculation shows that 
(4-87) lim zP^z) = 

z-^oo ZU 

oo 

Therefore, substitution of the expressions for / and into (I4.86P gives 



(4.88) h' 1 = - V e- 2t/3arcsinAo (l + 0(l/n 1 - 2 l^l)). 

We shall now calculate the recurrence coefficients in ( 11.101) and prove Theorem 11.21 Mul- 
tiplying the recurrence by p n (x)w(x) and integrating, we obtain by orthogonality 

/OO PfJ,() POO 

p 2 n (x)xw(x)dx = / f n (x)xe- x2 e i ^dx+ / p 2 n (x)xe- x2 e~ inf3 dx. 
-OO J —OO J flQ 

Integrating each of these integrals by parts, combining the resulting integrals together, and 
using the orthogonality of p n (x) and p' n (x), we obtain 

(4.90) A n = -h- l fMe-ti sinh(z7r/3). 
In terms of the RHP, 

(4.91) Pn{\vV2~n) = r u (A V2^) = (2n) n / 2 U ll (X ). 
In the region z(I) of the z-plane, 

(4.92) U n (X Q ) = S u {X )e n9M + S l2 {\ )e n9 - {Xo) . 

The main term is obtained here by substituting P\ q (\q) = lim 2 ^A P\ ( z ) f° r Here z tends 
to Ao along a path in z(I). Expanding P\ (z) we obtain that the main term of £/n(Ao) equals 

lim i ((a + a^V^-^- ty((3, + ^(1 - (3, e~™ Qf(/3)] e n ^ 

(4.93) W 

-i(a - a" 1 )^^ W-A e—C) + V>(1 + ft C)/(-0] e ^°)+ 2 ^j . 
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Recall that z — > Ao corresponds to £ — > 0. In this limit the ratio T>{z)jCf is given by 
(I4.67p . and the behavior of ip is determined by (14.701) . From the latter formula, as z — > Ao, 

(4.94) *l>((3, C) + Vi(l - /3, e-™()f((3) = e"^r(l - 0) + 0((z - A ) In(z - A )), 

(4.95) ^(-(3, e~**Q + + (3, C)/(-/3) = e"^r(l + /3) + 0((z - A ) \n(z - A )). 
A simple computation gives 



(4.96) <7 + (A ) = Aq — - - In 2 + 0+(A o ), 0+(A o ) = ztt \ ^(s)ds. 

1 JXo 
Substituting these equations into (14.931) yields 

Y n (\ V2ri) = (2n) n ^U u (\ ) = -(2n)"/ 2 P 00 2- n - 1 e n(A ^ 1/2) /5 

(4.97) x [(a + (A ) + a + (Ao)- 1 )F« + *(a + (A ) - a + (A )- 1 )F(-)](l + OCl/n 1 " 2 ^)) , 

= (8n) ±/3 (l - A 2 ) ±3/3/2 r(T/?)e ±n ^ (Ao) , P ra = e i/3 arcsin A ° . 

By (14.901) . (14.911) . this proves the formulas f II . 1 3 j) for A n in Theorem 11.11 The particular 
cases (11.151) . (I1.17P of the theorem follow from the simple identities: 

|r(±z 7 )| 2 = — \ , ar g r(-z 7 ) = - argr(^ 7 ). 7 G R \ {0}. 

7 smh 7T7 

In order to derive the formulas for B n in Theorem 11.11 and to prove Theorem 11.21 using 
Proposition ^. ll in the next section, we shall need more precise asymptotics for the coefficients 
x n , and also those for the coefficients /3 n , and 7„ of the polynomial p n {z). As usual, we 
compute them investigating the limit z —>■ oo of Y(z). By (12.31) . 

(4.98) ^-Hm.gg, M^^^ + ^t... 



z—*oo 



As 2 — > oo, we need to know asymptotics of Y(z) in the area outside the lenses (denote it 
A), which are given by the expressions: 

Y(zVZn) = (2n) na * /2 U(z), 
(4.99) U{z) = e nlas/2 (I + R^z) + 0(l/n 2 - 2 ^%P^(z)e< 9{z) - l/2)a \ 

z E A n Uoc. 

Set 



4 ±} = (8^) ±2/3 (l - Ag)^/(±/3)e^<*»\ L n = IH - L« 

(4100) = (Ao + V 1 "^) " ( A ° " V 1_A °) L " 

iV n = (a - i^/l-Ag) I« - (a + i^l-AgJ 

Let us compute R\(z) using (14.851) . Consider first the neighborhood U\ Q . Substituting 
Ai(x) given by (14.711) into (14.851) and calculating residues at a simple pole x = Ao, we obtain 
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the contribution to R\ from the neighborhood U\ : 



(4.101) 



27li /flrr, x — 2 z \ z I 2tx% 



n(l - A 2 ) \N n - 2if3 2 -2X f3 2 + iL n 



oo ' 



as z — > oo, where L n , M n , and iV n are defined in (I4.100p . 

To compute the contribution from the neighborhood Ui, we repeat the calculation now 
using Ai(z) from (14.761) . An additional complication is that the pole at x — 1 is of second 
order. We obtain, using (14.181) . 

1 2m 




/ 1 + 8p 2 ^ 4i/3 - 8ip 2 j^ + 8/5 • 



(4.102) x^l ' i_A0 , ' 1_Ao V !- A o I p 

-Ao 



-1 - 8/3 2 i±fa 



^3 

OO 



Here the dependence on {3 comes from the expansion (14.181) . 
Using (I4.19p . a similar calculation for C/_i gives 

>(-l) _ 1 f _ 1 f-i 1 , n/-«-2^ 



i? p = _ / = _ 1 _ _ + (z 

2m JdU-i x - z z \ z 

V°J ( "I " 8/3 2 ^ 4z/3 - 8*/3 2 ±±^ - 8/3^ 



l+Ao 
Ao 



(4.103) x-^ i_Ao , i_A ° V v 



^3 
OO 



3 x 64n 1 

Summing up the contributions (14.1011) . (14.1021) . and (14. 103ft . we obtain 

(4. 104) R 1 = R^ + + ^S Ao) • 

Substituting this into (14.991) . and using (14.161) and the expansions for z — > oo, 

(4.105) a(z) = 1 - 1 + -L + 0(z-*), g(z) = \nz-^- 2 + 0(z~ 4 ), 
we finally obtain from (14.981) 

(4.106) K n _ x = — — — <1 + — -z roT + TT-l v2+° 



0r(n-l)! I 4n(l-Ag) 2n 1 - A 2 



3A /3 2 - *L„/2 



2-2|J?/3| 



(4.107) /?„ = v / 2n^/3,/l-A 2 + ^- + O 



4n(l - A 2 ) \n 2 - 2 \W 
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In = n f - (3 2 (l - Aq) + i(3\Jl -XI+ , n * (3A§/3 2 + iN n /2 
(4.108) 1 ^ ( o) 

+2z/3^/l - A§[3A /3 2 - iL n /2] J + O 

where L n and iV„ are defined in (I4.100j) . The error terms here and in (14.97!) are uniform for (3 
in a bounded set provided only — 1/2 < 3?/5 < 1/2. Indeed, these terms in (14. 106j) - (14.1081) 
are those in (14.991) . at worst multiplied by a polynomial in (3 (independent of n) coming from 
the expansion (I4.16P of D(z). Now the error term in (14.991) is that from (14.841) which has the 
above uniformity property. A similar analysis holds for the error term in (14.971) . 

The combination of (I4.106P and (II. lip gives the asymptotics (I1.14p of Theorem 1 1.2 1 (see the 
remark following (15. 3p below), which reduce to (I1.16p . (11.181) in the particular cases (3 = i'j, 
and j3 = ij, A = 0. Thus, Theorem 11.21 is proved. 

We constructed a solution to the Riemann-Hilbert problem of Section 2 for n > uq, (3 in 
any bounded set of the strip —1/2 < 3?/3 < 1/2. By uniqueness, it gives the orthogonal 
polynomials via (I2.3p . On the other hand, the determinantal representation for the orthogo- 
nal polynomials shows that R(z) is an analytic function of (3. Furthermore, R\(z) is also an 
analytic function of (3 by construction. Thus, the error term in (14.841) is both analytic and 
uniform in (3 in a bounded set of the strip — 1/2 < 9ft/3 < 1/2. Therefore, it is differentiable 
in j3 (the derivative being uniform in (3). Hence, we easily conclude that the error terms 
in (I4.97p . (14. 1061) (14. 1081) have the same differentiability property. Alternatively, we could 
have deduced the differentiability of the error terms by noticing first that the asymptotic 
expansions of hypergeometric functions we used are differentiable in (3. 




5. Evaluation of the Hankel determinant 

Let us obtain the asymptotic expression for the r.h.s. of the differential identity (13. 8p . 
The expression for Yu(\o\f2n) was found in the previous section and is given by (14.971) . A 
similar representation holds for y 2 i(AoV / 2n): 

, N Y 2l (\ Q V2n) = (2n)- n / 2 U 21 (X ) = -(2n)- n l 2 ^-T^e^ 1 ^ (3 

(5.1) J^oo 

x [i(a + (\ ) - a + (A )- 1 )F(+) - (a + (A ) + a + (A )- 1 )F(-)](l + 0(l/n 1 - 2 ^)). 
Since our asymptotic formulas are differentiable in (3, we obtain from ( 14.971) and ( 15. II) : 



V' 



r ll y 21,/3 ~ r il,/3 r 21 - r H r i 
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(6.2) 



{^»V<-» - F«F<">' + O (^) } = -^YnY, 

-^4 2in(8n(i - A ° )3/v ( in ^)/°(^ 

As follows from ( 14.881) . 

[ln(x n x n ,)]^-2^ + o(-^). 



GAUSSIAN WEIGHT WITH A JUMP 



25 



Therefore we have for the last part of ( 13. 8ft : 
1 

7T 



e ^ S i n (tt/3) [Y n (^0)^21(^0)^ - >li (^0)^21 (>o) - (In >^>^-i)JjliiO*o)*2iO*o)] 

(5.3) 



We now turn to the evaluation of the rest of (13.81) . Care is needed with the estimation of 
x n . To obtain the asymptotics of K n from (I4.106P we need first to replace n with n + 1 and 
second, to replace Aq with Xq^^j. After a straightforward calculation, we obtain 



-n(ln x n x n _0^ - ( -2=i ) +2 [7^ - AX J 

(5.4) 



2n In iA + a/ 1 — An - 2/3 + 2iA W 1 - Xi + O 



Inn 

n i-2\m 



The sum of Q and ([53D yields 



In D n (/3) = 2in arcsin \ — 2(3 + 2iX n\/ 1 — Aq 

(JjL) 




<5 ' 5) -^^-^ /2 )-(^)'y°(^Sm) 

Since the error term here is uniform in (3, we can integrate this identity. 

Let fl be a bounded subset of the strip — 1/4 < 9ft/3 < 1/4 (cf. Section 2). The identity 
(15.5p was derived for /3 ^ f2 (see Sections 2 and 3). Note that the number of points (3 in fl, if 
any, is finite. Indeed, the function x 2 = x%((3) is a ratio D k /D k+ i of two analytic functions 
of /3, which are not identically zero because they are known to be positive for 9ft/3 = 0. Let 
us rewrite the identity ( 15. 5ft in the form f'((3) = 0, where f((3) = D n {(3) exp(— L r(n,t)dt) 
and where r{n,(3) is the r.h.s. of (13.81) . Since expression (15. 5p for r(n, (3) holds uniformly and 
is continuous for (3 G fl provided n is larger than some n (fl), and D n ({3) and its derivative 
are continuous, the function f((3) is continuously differentiable for all n > no(fl). Hence, 
f'(j3f) = for all Peto and n > n (fl). Taking into account that f{(3) = f(0) = D n (0) ^ 0, 
we conclude that D n (j3) is nonzero, and that the identity (15.51) is, in fact, true for all (3 e fl 
if n is sufficiently large (larger than no(fl)). 

Now set (3 = 0, and integrate (15.51) from to some (3 over (3. We obtain 



In !>„(/?) - lnD n (0) = 2nf3\n (iX + yjl- A 2 J - /? 2 + 2if3X n^l-X 

(5.6) 

/3 2 ln(8n(l - Aq) 3/2 ) - jf* In + /31n ^ + O ( 



n l-4|«0| 



In view of ( 11.61) Theorem 11.11 is proved. 
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7. Appendix 

The confluent hypergeometric function ip(a,c;z) is defined as a unique solution of the 
confluent hypergeometric equation 

(7.1) zw" + (c - z)w' -aw = 

satisfying the asymptotic condition 



(7.2) ^(a,e;z)~^B V n)41 J;r C K 

n=0 

37T 37T 

2^00, — -<argz<y. 

Here the standard notation 

/ \ i / \ / 1 \ / i \ r(a + n) 

(a)o = 1, (a)„ = a(a + 1) ■ ■ ■ (a + n - 1) = , — , n>l 

r(a) 

is used. The function ip{a, c; z) admits the following integral representation [3]: 

(7.3) c; z) = — — / ^(l + tY-^e^dt, 

r(a) Jo 

(7.4) — 7r < a < 7T, — — + a < argz < — + a. 
In (17.31) the inequality 

(7.5) > 

is assumed, and the branches of the functions t a ~ l and (l + t) c_a_1 are defined on the t-plane 
cut along the semi-axis (— oo, 0] and fixed by the conditions, 

(7.6) — 7T < argt < tt, — n < arg(l + t) < ir. 

When a = 0, the integration in formula ( 17. 3ft is performed along the positive semi-axis, 
< t < oo, and the formula gives ip(a, c; z) for the values of &rgz between — | and |. The 
increase of a from to tt leads to the clock-wise rotation of the ray of integration and, 
simultaneously, to the counter-clock-wise rotation of the domain in the (universal covering 
of ) punctured z-plane. This process produces an analytic continuation of i[)(a, c; z) to the 
domain 

71 37T 

(7.7) -<arg2;< — 

and is shown in Figure HI 

When z is in the domain (17.71) . the contour of integration is the lower side of the cut 
(— oo,0]. Similarly, the decrease of a from to —it leads to the contour-clock-wise rotation 
of the ray of integration and, simultaneously, to the clock-wise rotation of the domain in the 
(universal covering of ) punctured z-plane. This process produces an analytic continuation 
of ip(a, c; z) to the domain 

(7.8) - — <vcgz <-- 
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t-plane 



z-plane 



-1 



a = 



(a) 



* <argz<^ 
2 2 



-1 



71 

a= — 
2 



(b) 



< arg z < n 



-1 



a = 7i 



(c) 



| < arg z < — 



3jt 
2 



Figure 4. Analytic continuation of ^(a, c; z) to the domain | < argz < 



3- 



and is shown in Figure [51 When z is in the domain (17. 8p . the contour of integration is the 
upper side of the cut (— oo, 0]. 

When a goes beyond the interval (— 7r, 7r), equation (17. 3p provides the analytical contin- 
uation of the function ip(a, c; z) to the whole universal covering, C \ {0}, of the punctured 
z-plane C \ {0}. 

As z passes from one sheet of the universal covering C \ {0} to another, a loop around 
the interval [—1,0] must be added to the contour of integration in (17.31) . In Figure [6l we 
illustrate the process of the analytic continuation of ip(a, c;z) to the sheet marked by the 
condition 



(7.9) 



3n 5n 
— < argz < —. 



The domain (17.91) corresponds to the choice a = 2tt in (17. 3p with the contour of integration 
as shown in Figure [6)1. Assume that, in addition to (17. 5p . the inequality 



(7.10) 



3fc > 



holds. The loop in Figure [6H can be deformed to the sides of the cut [—1,0] with the following 
arguments of t and 1 +t in the integrand in the r.h.s. of (17. 3p : 



(7.11) 



argt|[_i, ]_ = argrj|[_i i0 ]_ 



-7r, argrj|[o i00 ) = -2tt, 



arg(l + t)|[ 



-1,0]. 



0, arg(l+i)|r 



-1,0]. 



arg(l + *)| 



[0,oo) 



-2tt. 



Here we use the notation [—1,0]+ and [— 1,0]_ for the upper and lower sides of the cut 
[— 1, 0], respectively. It follows from (17.111) that the function ip(a, c; z) on the sheet (17.91) can 
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t-plane 



z-plane 



-1 



a = 



(a) 



71 71 

— < arg z < — 
2 2 



a = 



2 



-1 



(b) 



- ti< arg z < 



oc = -7t 



(C) 



f <argz<-f 



Figure 5. Analytic continuation of ip(a, c; z) to the domain — < argz < — §. 



t-plane 
-1 




-1 



0C= 71 



a = 



5tc 



(a) 



(b) 



z-plane 



71 3Jt 

— < arg z < 

2 6 2 



3» <argz<I* 
4 6 4 



-1 



3tc 



(c) 



71 < arg z < 271 



-1 



a= 2ti 



(d) 



3tc 



5tc 



< arg z < 
2 5 2 



Figure 6. Analytic continuation of if) (a, c; z) to the domain ^ < arg z < ^f. 
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be represented as 



if>(a, c; z) = -r—^ / It I " 1 1 1 + tf-^e-^dt 

r(a) Jo 

(7.12) -^-^ e -^-a) / |t| a - 1 |l + t| c - a - 1 e -^ 
r(a) 7-1 



p -2ina /-oo o r 



r(a) 7 2 b 2' 

where we assumed the convention that arg \a\ = in any expression of the form \a\ b . Observe 
now that in the principal sheet, i.e. when — | < argz < |, the function ip( a ^ c 'i z ) is given 
by the original formulae (17. 3p - (17. 6p where one has to set a = (see also Figure Hh) • This 
yields the equation: 

I /*°° 7T 7T 

(7.13) ^ a ,c;z) = —J^ ^{l + t^e^dt, --<argz<-. 

Combined, equations ( 17.121) and ( 17. 13ft imply the following relation between the values of 
the functions i/;(a, c; z) and ip(a, c; e 2m z): 

(7.14) ^(a,c-e 2wi z) =e- 2mc ^(a,c-z) + -—(l-e- 2wi(c ~ a) ) / \t\ a - l \l + t\ c - a - l e~ zt dt. 

r(a) 7-i 

Note that this equation holds on the whole universal covering C \ {0}, i.e. for all values 
of arg z. Indeed, the function ?/>(a, c; z) has already been defined as a function on C \ {0}. 
The function if) (a, c; e 2m z) is a composition of the function ip(a,c;z) and the conformal 
isomorphism of the universal covering C\{0} determined by the mapping z — *> e 2m z. Finally, 
the integral over the finite interval [—1, 0] is an entire function and hence an analytic function 
on the universal covering C \ {0} as well. 

The integral term in (I7.14p actually describes the first canonical solution of equation (17.11) 
at z — 0. More precisely, equation (17. ip possesses a unique normalized holomorphic at zero 
solution, 



(7.15) (P(a,c;z) = J2 



(a)n z n 
(c) n n\ 



n=Q 

The integral representation of 0(a, c; z) links it with the integral from (17.141) . Namely, 

(7.16) <f>(a,c;z) = — ^ f f-^l - tf^'^dt 

Tare— a 



Vic) 



o 





\t\ a - l \l + t\ c - a - l e- zt dt. 



r(a)T(c — a) J_ x 
Therefore, equation (17.141) can be interpreted as the relation 

27T? 

(7.17) ^(o, c; e 2m z) = e~ 2 ^(a, c; z) + e™ - r0(a, c; z) 

r(c)T(l + a — c) 

where we used the classical formula r(s)T(l — s) = it / shi7rs. 
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It is worth noticing that the function 0(a, c; z), as an entire function of the variable z, is 
invariant under the mapping z — > e 2m z: 

(7.18) (f)(a,c;e ±27Ti z) =(f)(a,c;z). 

Our next task is to find another relation between the ^-and 0-functions. To this end, let 

(7.19) ${a, c; z) = e"^ (c - a V(c - a, c; e~ in z)e z . 

The function ijj(a,c;z) is well-defined by this equation as an analytic function on C\ {0}. 
In fact, i/)(a, c; z) is the second canonical at infinity solution of the confluent hypergeometric 
equation (17. ip . Its behavior at infinity is given by the asymptotic series 

(7.20) $(a, c; z) ~ £ ^~ a )f' a )^ 

n\z n 

n=0 

71 5tt 
z ->• oo, --<argz<y. 



Let us establish the relation between the functions ij)(a, c; z), ip(a, c; z), and <f>(a, c; z 

~~2 ^ arg «, , 



We start with the integral representation for ijj(a,c;z) assuming that ^ < argz < 



The latter implies that 

71 -in 37F 

- < arge z < — , 

and hence we can use the definitions ( 17. 3ft - (17.61) with a = 7r in (17. 19ft . This leads to the 
equation: 

^ „-i7r(c-a) p-oo o r 

(7.21) ^(a, c;2 ) = - / ir^l + t^Wdte?, — < argz < — , 

T(c - a) Jo 2 2 

where it is assumed that (cf. Figure Hb) 

arg t\ (-00,0] = -7r, arg(l + *)!(_«, = -7r, arg(l + t) | (— i,o] = 0. 

Changing the variable, 1 + 1 — > — t, and taking into account the above conditions on the 
arguments, we can rewrite (17.211) as 

^ -in(c-a) roo o r 

(7.22) V(a, C ;^) = r / (-l-ty-^i-ty-^dt, ^ < argz < ™ 

r(c-a)y_! 2 2 

where 

/ i -f\c — d — 1 1 ii i +|c— a— 1 m(c— a— 1) 

, 723) C— 1 — *J l(-i,oo] = + e 

/ j.\a— ll _ Ulo-1 / — 1 1 _ Ula— 1_— i7r(a— 1) 

(-*) I (—1,0] - 1*1 , (-*) |[o,oo) - 1*1 e 
This, in turn, leads to the equation 
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which, with the help of ( 17.131) and ( 17.161) . can be re- written as 

^a,c ]Z )=e- 2 ^~ a ^<P(a,c;z) 

(7.25) ^ 
-e- M ™ f7 ^y^(a, c; e" 2 ^), ^ < arg^ < ^. 

Note that the restriction ^ < argz < ^ can now be dropped; indeed, all the functions 
involved are analytic functions on the whole universal covering C \ {0}. Equation (17.251) is 
the relation between the functions vp(a, c; z), ip{a, c; z), and 0(a, c; z) we were looking for. 
Equations (I7.25P and (17.191) yield the formula 

(7.26) <f>(a, c; z) = r e^Ma, c; e~ 2ni z) + lM e ^ (c - a V(c - a, c; e~™ ' z)e\ 

T(c — a) r(a) 

which, with the help of the mapping z —>■ e 2m z, can be also written as 

(7.27) <f){a, c; z) = J^L e - i7r >(a, c; z) + ^r\e^ c ~ a ^{c - a, c; e in z)e z . 

1 (c — a) 1 (a) 

Observe that the operation z — > e~ 2m z brings equation ( 17.171) to the form 

(7.28) ij(a, c; e~ 2m z) = e 2 ^(a, c; z) - e™ -. -Ma, c; z). 

1 (c)l (1 + a — c) 

Using this relation in (17.261) . we can exclude ip(a, c; e~ 2m z) from the latter and obtain a 
"companion" equation to (17.271) . i.e., 

(7.29) 0(a, c; z) = , e CT >(a, c; *) + ^- e -^ c - a) ^{c - a, c; e~ in z)e z . 

1 (c — a) 1 (a) 

Note again that all the five relations (TT26|) . (173711 . dTTTTD . (TT^Bl . and (17^9]) hold on the 
whole universal covering C \ {0}, i.e., for all values of argz. 

Each of the equations (17.271) and (I7.29P allows us to exclude the function 0(a, c; z) from 
(I7.17p . This leads to the following two representations of the function ip(a, c; e 2m z): 

2tti 

(7.30) ^(a, c; e 2 ™z) = e'^^a, c; z) + e~™ a ^(c - a, c; e^)e 2 

1 (a)l (1 + a — c) 

and 

i/>{a, c; e 27ri z) = (l + e~ 2 ™ c - e - 2i7rc+2i7ra )ip{a, c; z) 

( 7 - 31 ) 2m 

1 (a)L (1 + a — c) 

Remark 7.1. Restrictions (17.51) and (17.101) can be lifted. The functions if) (a, c; z) and 0(a, c; z)/r(c) 
are in fact entire functions of the complex parameters a and c. For the function ip(a, c; z) 
this can be seen in the usual way by replacing, cf . [3] , the integration in (17.31) along the ray 
by the integration along the loop around this ray. Namely, one can rewrite (17.31) as 

i rooe~ la 

(7.32) 1,(a, c; z) = - ^.. TV . / ^(l + tf-^e^dt 



(1 -e 2 ™)T(a 



(o,+) 
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a = 





571 

4 



Figure 7. Contours of integration for ip{a, c; z) with arbitrary a and c. 



(7.33) 



2tt 



e-*™r(l - a) 



(o,+) 



^(l + t) 



c— a— 1 — zt 



(it, 



(7.34) 



7T 7T 

— h a < arg 2 < — 
2 & 2 



a, 



where the symbol /(^\ stands for the integration along the loop around the original contour 



in (17. 3p . In Figure [71 the new contour is shown for several different values of the parameter 



a. 



Equations (17.321) - (17.341) define the analytic continuation of ip(a } c; z) into the whole a- 
and c-complex planes (note that the integer values of a are removable singularities of the 
right-hand sides in (I7.32p . (17.331) ). Simultaneously, the validity of the relations (I7.30p and 
(I7.3ip is extended to all complex a and c. 

Also in the usual way (see e.g. [31J, Chapter 4, Section 5), the representations (I7.32p - 
(I7.34p can be used to ensure that the asymptotics (I7.2p are uniform for a and c belonging to 
any compact set of the complex plane. 

In the case of the function <j)(a, c; z)/T(c), the relevant analytic continuation is achieved 
with the help of the replacement in (I7.16P of the integration along the interval [0, 1] by the 
integration over the double loop: Pochhammer's loop (see, e.g., [3]). Alternatively, one can 
just use the relation (I7.27P and already established analyticity of the ?/>-function. 

Remark 7.2. In the standard literature, see e.g. [3], there is a convention to use equation 
(I7.29P for > 0, and to use equation (17.271) for Qz < 0. In fact, both equations are usually 
combined into a single formula: 



T(c) 



(7.35) 



T(a) 



T(c — a) 

-in(c-a)e j ( _ 



~ip(a, c; z) 



ijj(c - a,c;e lne z)e z 



e = sign Qz. 



The reason for this is to ensure that when — tt < arg z < ir the both -0-functions in the right- 
hand side have their arguments within the interval (— 37r/2, 37r/2) and hence the canonical 
asymptotics (17.21) hold for both terms. It is worth noticing, however, that there is no prohi- 
bition for using (I7.35p . say with e = 1 but for 3z < 0. 
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It is also worth mentioning one more time that, similar to equations (17.301) and (17.311) . 
equations (17.271) . (17.291) . and (17.351) are valid for all complex values of the parameters a and 
c 7^ — n, n = 0, 1, 2, . . . 
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